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Abstract 

In this work we investigate a phase field model for damage processes in two-dimensional 
viscoelastic media with nonhomogeneous Neumann data describing external boundary forces. 

In the first part we establish global-in-time existence, uniqueness, a priori estimates and 
continuous dependence of strong solutions on the data. The main difficulty is caused by the 
irreversibility of the phase field variable, which results in a constrained PDE system. In the 
last part we consider an optimal control problem where a cost functional penalizes maximal 
deviations from prescribed damage profiles. The goal is to minimize the cost functional with 
respect to exterior forces acting on the boundary which play the role of the control variable 
in the considered model. To this end, we prove existence of minimizers and study a family 
of “local” approximations via adapted cost functionals. 

AMS Subject classifications: 35A01, 35A02, 35D35, 35M33, 35M87, 35Q74, 49J20, 74A45, 
74D10, 74F99, 74H20, 74H25, 74P99; 
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1 Introduction 

Damage phenomena in elastically deformable solids and their analytical studies have received 
a lot of attention in the mathematical literature, e.g., 0 UHl [HI dSl ESI ESI ESI ED ESI SO]- 
Especially models which employ a phase field approach and incorporate higher-order terms 
were the focus in some recent works. In that case, an internal variable indicates the degree 
of structural integrity and, depending on the material and the scaling, may be defined as the 
volume or the surface density of microvoids or microcracks, respectively, as pointed out in m- 
This approach has also been utilized for approximations of surface discontinuities occurring 
in the displacement field of fracture models and turned out to be very useful for numerical 
implementations (see [DElEllEZ!). 
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One of the main difhcnlty for a rigorons mathematical investigation of the nnderlying PDE 
systems is that the damage variable is forced to be monotonically decreasing in time (irre¬ 
versibility) and, if possible, bounded in the unit interval. This kind of non-smooth evolution 
had motivated different concepts of weak solutions and regularization techniques in the literature 
(cf., e.g., [U [23l El] )■ However, to the authors’ best knowledge, a global-in-time well-posedness 
result for strong solutions with inhomogeneous boundary data was left open. Together with 
sufficiently strong a priori estimates such a result could be exploited to study optimal control 
problems typically arising in engineering problems focused on resistance against damage and 
failure. The following model problems with boundary control illustrate some practical exam¬ 
ples: 

- Suppose that a workpiece is exposed to external forces during an experiment and that 
certain parameters related to those forces can be controlled. A control problem could be 
to choose optimal parameters in order to prevent further damage in the material. 

- Related to the hrst scenario we might be interested in calculating additional forces not to 
prevent but to redirect crack spreading to non-critical components of the structure and to 
avoid complete failure. 

- Another problem might be the determination of external forces in order to deliberately 
induce a damage progression. For instance, it might be desirable to separate certain parts 
of the workpiece in industrial processes. 

By now, to the authors’ best knowledge, the mathematical contributions addressing those 
and related problems are inspired by the pioneering work [S] and employ fracture models to 
control the energy release rate of a single crack in a quasi-stationary setting by optimal shape 
design techniques, fibers or applied forces (see im [221 im EH] for more details). The cracks are 
explicitely modeled by non-smooth domains with or without non-penetration conditions for the 
deformation. A main issue consists in determining optimal forces or inclusions in the solid in 
order to cease crack propagation or to release as much energy as possible. 

In this paper we would like to advance a different approach for such control problems by 
utilizing a phase field model for damage. The kind of model under consideration was motivated 
by Fremond and Ned jar in m and is stated below. Under certain structural assumptions 
we are able to investigate well-posedness of strong solutions and existence of optimal boundary 
controls for a coupled evolutionary system describing damage processes in viscoelastic materials 
in two spatial dimensions. This enables us to investigate optimal control problems where we 
aim to control the damage phase field variable via external boundary forces. The cost functional 
will measure maximal deviations from desired phase held prohles. 

In the hrst part of this paper we study existence and then, for constant viscosity D, well- 
posedness of the following PDF problem: 

For a given time interval (0, T) and reference configuration D with boundary T and outer unit 
normal v, find {u, x) such that 


Utt - div(C(x)e(n) B(x)e(jJt)) = i 

in Q X (0, T), 

(la) 

Xt - Ayt -Ax + ^ + ^C'{x)£{u) : e{u) + f'{x) = 0 

in Q X (0, T) 

(lb) 

with the subgradient 



e <9/(-oo,0](Xt) 

in Q X (0, T) 

(2) 
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and the initial-boundary conditions 


u(0) = u°, utiG) = x(0) = X° 

in n. 

(3a) 

(C(x)e(n) -k D(y)e(nt)) ■v = b 

on T X (0, T), 

(3b) 

V(x + Xt)-J^ = 0 

on r X (0, T). 

(3c) 


Equation (jlall describes the balance of forces in the workpiece according to the Kelvin-Voigt 
rheology. The displacement field is denoted by u, the external volume forces by i, the linearized 
strain tensor by e{u) = ^(Vu -|- (Vu)^) and the stress tensor by a = C(x)e(tt) -|- B(x)e(tti). 
The first summand of a contains the elastic contribution whereas the second summand models 
viscous effects. The coefficient C designates the fourth-order damage-dependent stiffness tensor 
and D the viscosity tensor. The second equation (llbp specifies the parabolic evolution law for 
the propagation of damage described by the variable y under the constraint m, where the 
subdifferential of the indicator function h(_oo,o] : K —t B U { 00 } is given by 

{0} if xt < 0, 

[0, 00 ) if Xt = 0, 

0 if Xi > 0. 

The Laplacians — Ay and —Axt model diffusive effects of y and xt and have a regularizing effect 
from the mathematical perspective. For a mechanical motivation of system ([I|)-(l3]) by means 
of balance laws and constitutive relations we refer to [ni [13 [E]. In comparison to certain 
phase field models of damage used in the literature, the higher-order viscosity —Ay^ is also 
incorporated in the damage law (llbl) (cf. [HE])- If will help us to perform the (global-in-time) 
Second a priori estimate in Lemma 12.91 (see also the remark after the proof of Theorem 12.lip 
and, from the modeling point of view, it originates from an additional gradient term of xt in 
the dissipation potential corresponding to system o, which is given by 

R{ut,Xt) = + ^|Vytp -k ^D(y)e(ut) : e(?xt) -k/(_oo,o](xt)) dx. (4) 

We would like to give the following interpretation for the subgradient constraint ([2]): 

By introducing the free energy F to system m as 

F{u,x) = + \^ix)£{u) : s{u) -k fix)) dx, (5) 

we may rewrite m as 

0 G d^^R{ut, Xt) + dy^F{u, x) or, equivalently, ^ = -xt + Ay* - d^F{u, y). 
with a complementarity formulation for ([3), i.e. 

xi < 0, e • Xt = 0, e > 0. 

We comment that for a mechanical interpretation it would be desirable to force the values of 
the phase field variable x to remain in the unit interval [0,1]. Principally, that can be achieved 
by incorporating an indicator function in the energy ([5]) (or another suitable potential) or by 
using, if available, a parabolic maximum principle in the sense that x(0) € [0; 1] S'-®- H 
implies x(f) S [0,1] a.e. in 14 for almost all times t G (0, T). However, the first option would 


dl{-oo,o]ixt) = < 
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lead to the occurrence of a second subgradient ip G <9/[o_i](x) in the damage equation (llbh which 
would possibly violate uniqueness of solutions (a behavior observed in [7] for certain double 
inclusions; see also [39l Remark 2.18]) whereas the second option leads to serious complicacies 
in establishing the maximum principle due to the regularizing term —Axt in (jlbp . If the tensors 
C(-) and D(-) are constant for non-positive values we can at least provide a pragmatical solution 
(see also Example 12.ip : 

To this end, let {u, x) be a solution to the above problem and let the initial data satisfy x^ € 
[0,1] in n. Then together with the irreversibility constraints xt G (~oo,0] we find x £ (“Oo, 1] 
a.e. in X (0, T). Let us consider the pair (u, x+) with the pointwise truncation x~^ ■= niax{x, 0} 
which thus satisfies x'^ € [0,1] a.e. in x (0, T). Now note that (n, still fulhlls the elasticity 
equation (fTal) in 12 x (0,T) when replacing x by x^ because C(x) = C(x''') and B(x) = B(x''') 
for all values x- The damage equation pbl) when replacing x by x"*" is satisfied in the subset 
{x > 0} C 12 X (0,T). In the complementary part {x < 0} ^ 12 x (0,T) the damage evolution 
for x^ ceases and stays at the minimum, i.e. x"*" = 0 and thus (in an a.e. sense) {x'^)t = 0. 

The second part of this paper is devoted to an optimal control problem. A cost functional 
J will measure the maximal deviation of the damage variable x from given prescribed damage 
profiles at the final time T and/or at all times in [0,T] (AQ,An,As > 0): 

b) := -^llx - X(9 IIl°°(Ox(o,t)) + ^llx(T) - Xr||L°°(nx(o,T)) + “^ll^llL2{rx(o,r);R")- (6) 

A minimizer (x, b) of under the constraint that x solves system ([I])-© for some displacement 
u and admissible (later specified) boundary data b indicates an evolution which approximates 
Xq and/or xt best in the sense of J. It is also possible to replace x by x'^ oii tbe right-hand 
side of ([6]) in order to account only for the non-negative values of x- 

In the following we summarize the main results of our paper: 

- In Theorem 12.111 we will prove existence of strong solutions for system ([I])-([3]) and for a 
so-called /3-approximation in two spatial dimensions. In the latter case we replace the 
subgradient ^ in (l2|) by a smooth approximation ^^{xt) with /3 > 0. On the one hand this 
enables us to perform the a priori estimates in Lemma 12.91 while, on the other hand, the 
/3-approximation might be helpful for further studies such as optimality systems, numerical 
implementations etc. We emphasize that the existence analysis constitutes the main part of 
this paper and strongly relies on the two-dimensional Ladyzhenskaya’s inequality originally 
devised for the 2D Navier-Stokes equations (see [25] and the calculation ([32]) L 

- Continuous dependence on the data (u°, v^,x^, b, £) and, in particular, uniqueness of strong 
solutions for system ([ID-(I3D are proven in Theorem l2.12l fsee also Corollary 12.13p under the 
assumption of constant viscosity D. We also establish a priori estimates for the solutions 
in Corollary 12.141 These results allow us to define the solution operator and constitutes 
the fundament for the considered optimal control problem. 

- Theorem 13.61 reveals existence to an optimal control problem where the cost functional 

penalizes deviations of the damage variable from given damage profiles in the L°°-norm 
(see (l6|)). The strong solutions of system ([I|)-(l3]) will be controlled via external boundary 
forces. We prove existence of optimal controls by using the /3-approximation in the proof 
of Theorem 12.111 to define a family of optimal control problems. The minimizers or the 
optimal controls of the family of /3-approximating control problems converge in a limit 
process (along a subsequence as /3 0) to an optimal control of the original control problem. 
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In other words, we show that optimal controls for the family of /^-approximating control 
problems are for some /3 > 0 likely to be “close” to optimal controls for the original control 
problem. It is natural to ask if the reverse holds, i.e., whether every optimal control for 
the original control problem can be approximated by a sequence of optimal controls of 
the /3-approximating control problems. Unfortunately, we will not be able to prove such a 
“global result” that applies to all optimal controls for the original control problem. The 
reason for that lies on the non-convexity of the optimal control problems (both the original 
one and the /^-approximating control problems) and consequently on the non-uniqueness of 
the optimal controls. However, a “local” result can be established by introducing so-called 
adapted optimal control problems in Theorem l3.ini 

Let us recall some already established results in the mathematical literature of phase field models 

for damage/gradient-of-damage models: 

- Local-in-time well-posedness of strong solutions for damage-elasticity systems with scalar¬ 
valued displacements, homogeneous Dirichlet conditions for the displacements and linear 
dependence of C and B on y is proven in HE] and in mm for one-dimensional models. 
Since a degenerating elastic energy with respect to y is considered and enhanced estimates 
are established for small times, the solutions are obtained locally in time. In the work [3| 
such rate-dependent damage-elasticity system with the additional second order term xtt 
has been coupled with an equation for heat conduction incorporating highly non-linear 
terms. Existence, uniqueness and regularity results are established locally in time. 

~ Rate-independent gradient-of-damage models are explored in m and in subsequent pa¬ 
pers, e.g., [301 132] • In that case the corresponding dissipation potential R is positive, 
convex, and positively-homogeneous of degree 1. From the modeling point of view it 
means that the damage progression is considered on a faster timescale than the acting 
of the external forces. The authors considered non-smooth domains and employed weak 
notions referred to as energetic formulation in order to prove existence of solutions. The 
degenerating case where the material may loose all its elastic properties due to heavy 
damage is also studied. Further cases involving nonlinear r-Laplacians with r > 1 or even 
r = 1 instead of the classical Laplacian in the damage equation are investigated in |33l 03] , 
where also higher temporal regularity is shown. 

- A weak notion for rate-dependent damage models coupled with Cahn-Hilliard equations 
was introduced in |18] for quasi-static balance of forces and in m with inertial effects and 
without the viscosity term in (|lal) . Existence of weak solutions is proven there for non¬ 
smooth domains and mixed-boundary conditions for the displacements whereas uniqueness 
is left open. 

- Well-posedness and vanishing viscosity results for damage models with (nonlocal) higher- 
order s-Laplacian are established in |23] (see also [23] for vanishing viscosity results for 
damage models with a regularized nonlinear g-Laplacian in non-smooth settings). The 
authors study both viscous (rate-dependent) and, via the vanishing viscosity limit, rate- 
independent PDF systems for damage. Existence of solutions is provided in both situations 
where in the latter case a novel energetic formulation making use of arc-length reparame¬ 
terization in order to describe the behavior of the system at jumps is utilized. In particular, 
the existence results cover the case of the standard Laplacian in two spatial dimensions. 
Furthermore, the uniqueness problem has been solved under special conditions, where in 
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the case of two or three spatial dimensions the s-Laplacian is assumed to be of higher-order 
than the classical Laplacian. 

- Coupled thermoviscoelastic and isothermal damage models incorporating p-Laplacian op¬ 
erators are analyzed in [39] (see also [iQ] for the full heat equation including all dissipative 
terms and |2r)| for damage-dependent heat expansion coefficients). In those works homo¬ 
geneous Dirichlet boundary conditions for the displacements are assumed. Uniqueness is 
shown in the isothermal case by adopting p > n and by dropping the irreversibility con¬ 
straint dSj). Existence results for the corresponding rate-independent thermoviscoelastic 
damage models are proven in the recent paper [26j . 

Structure of the paper 

Section [2] is devoted to the well-posedness problem of system m-m- We state the precise 
assumptions in Subsection 12.11 and introduce time-discretized and /3-regularized approximations 
of ([I])-([3|) in Subsection 12.21 The existence proofs are carried out in Subsection 12.31 firstly for the 
time-discretized and then, by a limit analysis, for the time-continuous versions. In the final part 
of that section, i.e. in Subsection 12.41 we prove continuous dependence on the initial-boundary 
data. Then, equipped with the well-posedness result, we state the announced optimal control 
problem in Section [SJ We prove existence of optimal controls via /3-regularization in Subsection 
13.11 and their approximation by means of an adapted cost functional in Subsection 13.21 

2 Analysis of the evolution inclusions 

The approach presented in this work combines two different approximation techniques to obtain 
existence of solutions for system m-m- semi-implicit time-discretization and regularization of 
the subgradient ^ in ([2]). At first we will tackle the existence problem for the time-discrete and 
regularized system in Lemma 12.71 By passing the discretization fineness to 0, solutions of a 
time-continuous regularized system are obtained in Theorem [2TT](i) . In the final step, a further 
limit passage leads to solutions of the desired limit system (see Theorem 12.111 (ii)). Then, we 
conclude this section in Theorem 12.121 with a uniqueness and continuous dependence result. 

2.1 Assumptions and notation 

Throughout this work, we adopt the following assumptions: 

(Al) n C M"- with n € {1, 2} is a bounded C^-domain. The boundary is denoted by T and the 
outer unit normal by v. 

(A2) The damage-dependent stiffness tensor satisfies C(') = c(-)C, where the coefficient function 
c is assumed to be of the form 

c = Cl -I- C 2 where Ci € (^^’^(M) is convex and C 2 G C'^’^(M) is concave. 

Moreover, we assume that c,c'^,C 2 are bounded and as well as 

c(x) > 0 for all x G M. 

The 4**^ order stiffness tensor C G £(M”y^; is assumed to be symmetric and positive 

definite, i.e. 

Cijik = Cjiik = Cikij and e : Ce > 7/|ep for all e G M”y^ (7) 


6 








with constant rj > 0. 

(A3) The damage-dependent viscosity tensor satisfies B(-) = d(-)D, where the coefficient func¬ 
tion d satisfies d € Moreover, we assume that d and d' are bounded and 

d(x) > r/ > 0 for all x G M and fixed rj > 0. (8) 

The order tensor D is given by D = fj,C, where /r > 0 is a constant. 

(A4) The damage-dependent potential function / is assumed to be in / G 

In the assumptions above should be understood as the space of differentiable functions 

on M whose derivatives are Lipschitz continuous. 

Example 2.1 Letc G 1]) be any given non-negative function with (c)'(O) = 0. Then 

there exists an extension ofc to the entire real line M (the extension is denoted by c) such that 
c(x) = c(0) for all X < 0 and assumption (A2) is satisfied for c and a given tensor C with Q. 

Proof of Example 12.11 We define the convex function ci and the concave function C 2 such 
that c = Cl -|-C 2 on the compact intervall [0,1] via 

Ci(x) := c(0)-I- J max{(c)"(T),0}dr^ ds, 

C 2 (x) := J min{(c)"(T),0}dr^ ds. 

To extent the functions Ci and C 2 from the domain [0,1] to M in accordance with (A2) we 
have to be careful because the extensions should have bounded derivatives as well as bounded 
second derivatives whereas the sum of the extensions should also be bounded. We provide the 
following construction: 

Since (ci)' > 0, (ci)" > 0 as well as ( 02 )' < 0 and ( 02 )" < 0 on [0,1] it hold 

(ci)'(l) = max |(ci)'(y)| =: Ai, -(c 2 )'(l) = max |(c 2 )'(?/)| =: A 2 . (9) 

ye [0,1] ye [0,1] 


In the case Ai < A 2 we may extend Ci and C 2 to M as follows (for readers’ convenience we keep 
the integrals explicitly): 


ci(x) 


C2(a;) 


ci(0) 

^ ci(x) 

fi -^1 + (^2 — Ai)^-j^ ds 

ci(l) + /i Ai -|- (A 2 — Ai)^-j^ ds -|- A2 (x — (1 -|- (5)) 

02(0) if X < 0 , 

< C 2 (x) if 0 < X < 1 , 

^02(1) — A2(x - 1) if 1 < X. 


if X < 0, 
if 0 < X < 1, 
if 1 < X < 1 -|- (J, 
if 1 -|- J < X, 


Here, d > 0 can be chosen arbitrarily. Due to ([9]) we observe that Ci G C^’^(M) and that Ci is 
convex in the entire real line M. Furthermore, C 2 G C^’^(M) is concave in M. The particularity 
of this construction is that Ci and C 2 have only linear growth for large x > 0 with derivatives A 2 
and — A 2 , respectively. Thus the sum ci -|- C 2 is bounded. 

Hence we observe that c'^,C 2 and Ci -|- C 2 are bounded and that the properties in (A2) are 
fulfilled for c, Ci and C 2 . The case Ai > A 2 can be treated analogously. □ 
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Remark 2.2 (i) The non-degeneracy condition dS]) prevents the material from complete dam¬ 

age, i.e., even the maximal damaged parts (the region with x < Oj exhibit small viscous 
properties. 

(a) The assumption D = piC in (A3) is needed in the proof of Lemma \2.1\ in step 2 in order 
to perform a regularity argument based on a transformation. It has already been employed 
in the mathematical literature (see \^ ). 

For later use, we define the solution space U x X, where U denotes the space of the displace¬ 
ments and X the space of the damage evolutions given by 

U := (10a) 

X := H^{0,T;H^{n)). (10b) 

The space of boundary controls B is defined as 

B := L^{0,T; n H^{0,T; 

We also introduce the sets for brevity 

Q:=nx{0,T), S:=rx( 0 ,T). 

Finally, let us mention that we make frequently use of the standard Young’s inequality 

ab < 5af + —h^ for all a, 6 € M and all <5 > 0 
4o 

where <5 > 0 will be chosen when necessary and we write Cs ■= Moreover, the symbols C, 
C, D, rj and 6 will denote positive constants throughout this work. 

2.2 Notions of solution 

Let us consider two approximations of system dU-Q: a regularized version where the indi¬ 
cator function /(_oo,o] ™ d21) is replaced by a suitable smooth function Ip, (3 G (0,1), and a 
time-discretized version of the regularized system. To this end, we introduce the following 
regularization: 

Definition 2.3 (/3-regularization) Let the family of functions {L/ 3 }/ 3 g(o,i) ^ (7^0(1^) denote a 
regularization of the indicator function /(_oo,o] following sense: 

(i) Ipi < 1/32 pointwise in M for every /3i,/32 € (0,1) with /3i > (32, 

(a) Ipfoo pointwise in [0, oo) as /3 fO, 

(Hi) Ip{x) = 0 for all X < 0 and all (3 G (0,1), 

(iv) I'p{x) > 0 for a.e. x G M and all f3 G (0,1). 

We may also write := Ip in the following. 


Remark 2.4 In particular, we may choose the Moreau-Yosida approximation given by (see 
Lemma 5.17]) 


1/3 (x) = inf 

ygR 



if X <0, 
if X > 0. 


Let us mention that also -approximations may he chosen for {Ij^} especially in view of opti¬ 
mality systems for optimal control problems (see Chapter 5]). 


Definition 2.5 (Strong solutions) For system ©-© and their approximations we introduce 
the following notion of solutions: 

(i) Time-continuous limit system (t = 0, 13 = 0). 

Let the data ,v^ ,x^,b,£) be given. A solution of the time-continuous limit system is a 
pair of functions {u,x) & Id x iY safisfying ©-© in an a.e. sense and for a subgradient 

(ii) Time-continuous (3-regularized system (r = 0,(3 > 0). 

Let the data {u^,v^,x^,b,£) be given. A solution of the time-continuous (3-regularized 
system is a pair of functions {u,x) ^ Id x X with u{0) = u^, dtu{0) = and x(0) = x^ 
such that 


utt - div(C(x)£(^x) +D(x)e(nt)) = £ a.e. in Q, 

Xt - Axt - Ax + ?/3(xt) + ^C'{x)£{u) ■ £iu) + fix) = 0 a.e. in Q, 
(C(x)e(u) + D(x)e(ttt)) • a.e. onY, 

V(x + Xz)'Z^ = 0 a.e. on S. 


(lla) 

(llb) 

(llc) 
(lid) 


(Hi) Time-discrete (3-regularized system (t,/3 > 0). 

Let {0 ,t,2t, ... ,T} denote an equidistant partition of [0,T] with discretization fineness 
T := T/M and M € N. Furthermore, let the data iu^,u~^,x^), {b^}k=o,...,M as well as 
{£^}k=o,...,M be given. A solution of the time-discrete (3-regularized system is a sequence 
{n^, X^}fc=o,...,M of functions G R^(11;R”') and x^ G such that 


_ 2^fc-i + ^k-2 


.yk y,k 1 y^k y,k l 

A, A. A A 


- div(c(x'')e(n'^)+B(x")£( 

Ax'' + ?/3( 




= & a.e. inLl, (12a) 




T 

+ ^(cUx"') + fix^-f)Csiu^-f : e(n^-i) + /'(x") = 0 


)eiuf + I]>ixfE 


u — u 


k-l 


V = h^ 


.yk - yk 1 

V(x" + ^—^ 


V = 0 


a.e. in Ll, (12b) 
a.e. on F, (12c) 
a.e. on T (12d) 


for all k = 1,... ,M, where c = Ci + C 2 denotes the convex-concave decomposition from 
(A2). 
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Remark 2.6 If we assume Vx^ ■ v = 0 a.e. onV we even obtain 

Vx ■ n = Wxt ■ n = Q a.e. on S 


( 13 ) 


instead of pel) or (llldl) and for all k = 1,... ,M 

_ .^k —1 

= ^-z/ = 0 a.e. on r (14) 

T 

instead of ()12dl) . 


2.3 Existence of solutions 

2.3.1 Existence for the time-discrete regularized system 

At first we are going to show existence of time-discrete solution according to Definition 12.51 (hi). 
Let r > 0 and /3 > 0. To enhance readability, we will mostly omit the subscripts r and /3 in ^ 
and Xr, 0 - 

Lemma 2.7 Let the data G x° G G R^/2(r;M”) and G 

L^(D;M”') for k = 0,..., M be given. Then, there exists a strong solution {u^,x^}k=o,...,M 
of the time-discrete system in the sense of Definition \2. 5\ (in). 


Proof. Starting from the initial values (uP,u~^,x^) with u~^ := — tv^ we are going to 

construct {n^, x^}fc=o,...,M by a recursive procedure. To this end, we decouple the discrete PDE 
problem into two distinct elliptic problems such that x^ is obtained from x^~^ u^~^, while 

is gained from u^~^, x^^ b^ and . 

Step 1: establishing equations (jl2bll and (I12dp 

Let us define the functional F : H^{Q) —M by 


nx) := 


^ (^|Vxl^ + ^Ci(x)Ce(n^ ^) : e{u’‘ ^) + ^^(x^ ^)xCe{u^ ^) : s{u'^ ^)) dx 

+ [ [fix)+ Tlg{{x-X^~^) t-~^)) <ix + ^ f ix-X^~^)r~^ 

J J 


^27, 


V(x-x''-')r-' 


dx 


By the direct method in the calculus of variations, we obtain the existence of a minimizer of F, 
which will be denoted by x^- 

The Euler-Lagrange equation for the minimizer yields (I12bl) in a weak form. By noticing 
that ()12bp is a elliptic equation for x^ with right hand side in L^(D), we conclude x^ G H‘^{Xl) by 
elliptic regularity results for Neumann problems (see, e.g. m Theorem 2.4.2.7] and remember 
that r is a C^-boundary by Assumption (Al)). 

Step 2: establishing equations (ll2aD and (|12c|) 

Given the functions x^ G R^(D), G b’^ G ^ 2>2(D;M"), 

we obtain a unique weak solution G R^(D;]R”) of the linear elliptic system (ll2aD via the 
well-known Lax-Milgram theorem (remember the assumption D = fiC from (A3)): 


((^^c(x'')+ rd(x'')/i)Ce(a) : e(C) + u • c) dx 
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- rdiv(d(x*^)De(7x'=-i)) + 2u^-^ - • Cdx + j 6^ • Cdx (15) 


holding for all ( € 

Now we use a modification of the regularity argument in |2Ul Proof of Lemma 4.1] and make 
use of the C^-regularity of P (see (Al)): 

If we consider the test-function ( = (r^c(x^) -|- rd(x^)l^) where ip € iL^(ll;R") is another 
test-function (remember that c > 0 and d > 77 > 0, see (A2)-(A3)) the linear elliptic system (|15p 
rewrites as 


a{u’^,‘p) = {q.,<p)m, ip e 


(16) 


with the bilinear form 


a{u,ip) := / Ce{u) : £{ip) dx 


and the right hand side q G given by 


(q, 

where R and N are defined as 


/ R- if dx + 
Jn 



__ r^c^(x^) + r/rd^(x^) ^ fc.p / fcx - rdiv (d(x^)D£(n^~^)) - + 2u^-^ - 

T^c(x^)+ rd(x^)/U ^ t2c(x^)-F rd(x^)/x 

bk 

N ;= -. 

T^c(x^) + Td{x^)p 


Note that N G ^ ^ LP{n-,R'^) for all p G (1,2), since e(u^) G L2(11;M^x^) and 

Vx*^ G LP{n-R^) for all q G [l,+oo). 

In particular, q G for all s G (1,3/2). We gain G Lf*(ll;M"') by apply¬ 

ing the lower Sobolev -regularity result from [HI Theorem 3.4.5 (ii)]. This, in turn, implies 
£{u^) G L^*(r2;R”) with the fractional critical exponent given in this case by 2* = > 2 

(see, e.g., [35l Theorem 6.7]). We obtain R G L^(n;R"). The i7^-regularity result [H Theorem 
3.4.1] applied to the linear elliptic system (11611 shows G 77^(11; R”). Thus (jl2ap is shown. □ 


2.3.2 Existence result for the time-continuous system 

The aim of this section is to provide existence of strong solutions in the sense of Definition [23] (i) 
and (ii). To this end, several a priori estimates for the time-discrete solutions will be established. 
The estimates will be used for the time-continuous limit analysis and for the optimal control 
problem in Section [3| 

We assume that the initial data (u^,v^,x^) satisfy 
v° G 77^(D;R^), 

X° G G H^(n) lVviy = 0 a.e. on T} 
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(17a) 

(17b) 

(17c) 







and the external forces {b, t) are assumed to be in the following spaces: 

^ € ^^(0,T;M")). (18) 


For the moment, let us consider some approximations 

{6A}Ae(0,l) CCl’l(0,r;/7V2(r.Kn)^ 

{4}Ae(0,i) CC'°’Ho,r;L2(b!;M-)) 

of the the initial velocity and the external forces b and I such that (e.g. construction via 
convolution) 


^ strongly in M”), (19a) 

bx^b strongly in L2(0,r;F^/2^r;R”)) nF^(0,r;L2(r;R”)), (19b) 

strongly in L^(0,T;L^(n;R"')) (19c) 

as A 4- 0. Let us define the time-discretizations b^^ and ^ by 

b’l^^-.= bx{Tk\ £’:^x--=hirk). 

For a sequence {h^}k=o,...,M where G {u^ x^^^ >} ^ dehne the piecewise constant 

and linear interpolation as 

h{t) := h{t) := , h{t) := 


h{t) := 


t — {k — l)r 
r 


h^ + 


T 


> for t € {{k — 1)t, kr], (20) 


The left-continuous and right-continuous piecewise constant interpolation for a given time point 
t is denoted by 


tr ■= rk for T{k — 1) < t < rk, 

:= rk for rk <t < T{k + 1). 

For notational convenience, we define the time-discrete velocity field and their interpolations by 

'^r,p '■= ^ for A: = 0,..., M and Vt,p,Vt,p, Vr,p by ([20]). (21) 

As a hrst result, we prove convergence of the discretizations of the given data. 

Lemma 2.8 There exist subsequences r*, | 0 and ), 0 as A: t oo such that 
—7> strongly in 

Km ^ b strongly in L\0,T-, n H\0,T-, L\T-R^)), 

Km ^ strongly in L^(0, T; L^(n; R"')) 

as A: t oo. For readers’ convenience we set br^ := brf.^Xk '■= ^TkM omit the subscript 

k. Then the statement above reads as ^ , b^ ^ b and £ as r ), 0. 
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Proof. For every fixed A > 0, we find 


(22a) 

(22b) 


br,x bx strongly in L‘^{0,T-, n 

^r,A —>■ strongly in L^(0, T; L^(fl; M"")) 


as r 4, 0. Indeed, the first convergence in ()22ap follows by exploiting the Lipschitz continuity 
of bx € C'°’^(0,r;Ff^/^(r;M”)). Property ()22bh can be proven with a similar argument. The 
convergence b^^x bx in the Ff^(0, T; L^(r; M"'))-norm follows by the fundamental theorem 
of calculus for AC-valued functions where X := L'^(T'W^) and by the Lipschitz continuity of 
dtbxGC^’Ho,T-,xy. 



dtbx{t) 




bxjU) - bxiU) 

T 


dtbx{t) 


2 

X 


dt 



{dtbx{s) - dtbx{t)) ds 


2 

X 


dt 


^ 0 J 


|s-t|ds)^dt 

< CTr^. 


The claim follows by using the convergences (|19D . (1221) and the following “e/2”-argument: 

For a given e > 0 we may choose a small A > 0 such that ||6 — 6 a|| < e/2 due to ([T9]). For 
every such A we choose a small r > 0 such that ||6 a — ^r,A|| < e/2 due to ([221) . In consequence 
we find for every e > 0 small values A > 0 and r > 0 such that ||6 — 6r,A|| < e. □ 


Lemma 2.9 (A priori estimates for the time-discrete system) The following a priori es¬ 
timates hold for strong solutions of the time-discrete system given in Definition \2.5\ (Hi) (recall 
that m implies = dtUr^p ): 

(i) First a priori estimate: 

There exists a constant C > 0 which continuously depends on 

C = Cdln^llj^l, ||u°||j;^ 2 , llx^ll/fl, ||^||L 2 ( 0 ,r;L 2 (r;R"))i ll^llL 2 ( 0 ,r;L 2 )) (23) 

such that for all t, (3 > 0 

ll''^T,^l|j^l(0,T;iTl)nVKl.°°(0,T;L2) < C, llXr,^ IliTpO.T;//!) < C, 

IIiIt,/3|Il°°(0,T;H1) < C, IIX.^^^IIl°°(0,T;H1) < C, 

lkT,/3||L°°(0,T;//l) < C, IIXt,/3|Il°°(0,T;H1) < C. 

(a) Second a priori estimate: 

There exists a constant D > 0 which continuously depends on 

D = iA(||tt°||j^2, ||x°||iy2, ||^||2,2('o^'r;//l/2(r;]R"))nHl(0,T;L2(r;R"))) PllL2(0,r;L2)) 
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such that for all r, f3 > 0 

lk'r,/3||//l(0,T;_H'2)nWl’°°(0,T;//l) < D, 

ll^r,/3||L2(o,T;H'2)nL°°(0,T;//l)nHl(0,T;L2) < D, 


IIXT,/3||//l(0,r;_H'2) < D, 
IIX^^^IIl°°(0,T;H2) < D, 

IIXr,/3llL°°(0,T;fl'2) < D, 

ll'^/3(^tXT,/3)||L2(o,r;L2) < D. 


Proof. We will omit the subscript r and /3 in the time-discrete solutions. 

To (i): In the following, we make use of a combined convex-concave estimate for: A convexity 
estimate for ci and concavity estimate for C 2 yield: 


Ci(x"-')-Ci(x'=)>c;(x")(x"-i-X^) 

C2{X^-^) - C2(x") > c'2(x''-')(x''-' - x'^)- 


Adding them shows 

c(x"-') - c(x") = (ci(x''-') - ci(x'=)) + (C 2 (x''-') - C 2 (x'=)) 

>(c'i(x'=) + c'2(x"-'))(x'=-'-x'=). 

By using this combined estimate and the positivity of C, it holds 
c(x'=)Ce(u'=):e(n'=-n'^-i) 

= ^c{x’')(Ce{u^) : e(u^) - Ce{u’^-^) : e{u^-^) + Ce{u^ - u^-^) : e{u^ - 
> £(|!lce(u^) : £{u^) - : e(u^-^) 

+ i(c(x"-i)-c(x"))C6(u^-^):e(u^-i) 

+ i(c;(x") + c'2(x'=-'))(x"-' - x'^)C£{u^-^) : e(u"-i). (24) 


Now, by testing equation (I12al) with integrating over Q, summing over the time 

index k = 1,... ,tr/T, integrating by parts and using (|12cp . we obtain (remember that 

yk _ (^yk _ yk-ly^-^ 


*r/r U/t 

^ / (/-/-^)u^dx + ^ / c(x'=)Ce(n^) :e(u^-n^-i)dx 

k=l k=l 


+ / / d(x^)De(u) : e(u) dxds 


0 Jn 


0 Jn 


l-vdxdiS+ / / b-vdxds. 


0 Jr 
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Applying elementary estimates including the convex-concave estimate (|24p . Korn’s and 
Young’s inequality and the trace theorem ^ L^(r;M"') yield {r],S,Cs > 0 are 

constants) 


\\\v{t)\\l2 - ^\\v°\\l2 + : e(u(t)) dx - ^ : e(u°) dx 

^-V-^ 

>0 by using (A2) 

+ [ [ x(ci(x)+C2(x))(-<9tx)Ce(u) : e(u)dxds-Fr/||e(u)||2 2^p-^.^2^ 

t/ 0 t/ 

^ ^sM\\l^{0,U;L^) + ^'5ll^llL2(0,KY2(r;R")) +'^II'^IIl2(0,Y;H1)- 


Testing equation (I12bp with ~ integrating over Q, summing over the time index 

k = 1,... ,tr/T, integrating by parts and using (|12dp . we obtain 


tr/r 




k=l ' 


UIt 

+ E 

k=l 
rtr 




{x^-X^-^)Ax + 


/o 


O (ci(x) + C2(x))9fXC£(u) : £{u) dxds 


+ 


[ [ f{x)dtxdxds 

Jo Jn 


= 0 


By using the monotonicity of (see Definition 12.311 . we get ^^ (x^ ~ X^ > 0- 

Together with elementary convexity estimates, the Lipschitz continuity of /' (see (A4)) 
and Young’s inequality, we find 

lia.xlli.(„,i„L.) + liva,xlli.,+ l||vx(t)iii. - l|ivx“ii'i, 

+ / / ^(4(x)+C2(x))5txCe(u) : e(u)dxds 

Jo Jn ^ 

— ^L'^{0Jt\L‘^) 

To proceed, we consider the calculation 


Adding dlY, 


o \j o — 

^ / / X^txdxds - ^||x°||i2 

Jo Jn ^ 

-,+C,||x|li.(oy,^^.) 


^Ilx°lli2. 


1^2 on both sides in 




and using (f27p on the right-hand side, we 


(27) 

find 


■ ll^tXlli2(oy,.^2) + l|Vatx|li2(oy,.^2) 


+ 


^llx(t)ll^i - hlx^Wm 
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+ [ [ ^(c'i(x) + C 2 ix))dtxCe{u) : e(:u)dxds 

Jo Jn ^ 

s + !)■ (28) 


Adding ([281) and (1251) . and choosing (i > 0 small, we see that the term 

[ ! \ (‘^Ux) + C 2 (x)) 5 txCe(:u) : e{u) dx ds 

Jo Jn ^ 

cancels out in the calculations and we obtain 


II \ /PI j_j 11 V- \ /II _n_ II \ ' 11 j_j \\j 

<(7^1 + + ||^^°||i2 + llx^ll^i 




h ds- 


Korn’s inequality yields 


We thus obtain from 


k(^) 11^2(0,> ^lbllL2(0,t^;_H-i) ll'*^llL2(0,i^;L2)- 


\Ht)\\h + llxWllii + Il^lli2(0j^.^i) + l|5iXlli2(oy^,^i) 

< c(l + Wu^Wjji + ||v°||^2 + Wx^Whi + Plli2(o,t^;L2) + \\^\L' 2 { 0 ,U;L^^■,R^))) 

+ C'^ (PIIl 2 + ||xp 2 ^ ds. 

We end up with the desired estimates in (i) by using the discrete version of Gronwall’s 
lemma and 


u 


2 

Hi 


«" + 


/ 


V ds 


Hi 


<ci(lh 


|2 _L 

Ihi + 


l^llL2(0,t^;Hl)) 


afterwards. 


To (ii) - local-in-time estimate; 

At first we are going to show the a priori estimates in (ii) for small time. In the next step 
global-in-time estimates will be derived. 

Testing equation (|12ap with —Tdiv^c(x^)C£(w^) + d{x^)Y)e( ^*‘~'!f integrating over 

n in space and summing over the time index k = 1,... we may write the result in 

the following way 

f f —dtv ■ div(c(x)Ce(’u) -|- d(x)De(T)) dxds 
Jo Jn 

' -V-' 

= :Ti 
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+ / [ ^1 div(c(x)Ce('u)+ d(x)D£(t;))|^dxds 

Jo Jn ^ 


+ 


10 Jn 


-ixWx • Ce{u) + c(x) div(Ce(u)) + d'(x)Vx • De(?;) + d(x) div(De(?;)) 


dx ds 


= :T 2 


10 Jn 


—i ■ div(c(x)C£(ti) + d(x)De(x)) dx ds . 


(30) 


<C'«IKII^2(i2)+'5|l 'iiv(c(x)C£(«)+d(x)De(i;))||2 2(^2) 


Note that the second summand and the third summand, i.e. r 2 , are identical. The splitting 
will simplify the calculations. 

Testing equation (|12bp with —A(x^ — X^~^), integrating over Q. in space and summing 
over /c = 1 ,... ,trlT^ we obtain 

- r [ SixASixdxds + IIAatxIli^.oi .^2) +V / Ax"(Ax'^-Ax"-')dx 
^ Jo Jn , V . T, ; 

^ ^ V_/ 

= l|Vaxlli2(,2) by using II14II 


'0 JQ. 


?/3(5tx)A5txdxds- 


/o Jn 


^(ci(x) +C 2 (x))Ce(n) : £(u)Aatxdxds 


= :T3 


=.Ta 


[ [ /'(x)A9ixdxds 

Jo Jn 


=.n 


= 0 . 


(31) 


In the following, we are going to estimate Ti,... ,Ti and conclude the claimed a priori 
estimates thereafter: 

- To (Ti): Integration by parts in space yields 


Ti = 


f f £{dtv) : c{x)C£{u) dx ds + f f e(cltu) : d(x)D£(t’) dx ds 
Jo Jn Jo Jn 


= :TT^ 


= :Tp^ 


- [ f dtv ■ ((c(x)Ce(?x) + d(x)De(u)) • u) dx ds 
Jo Jr 

'-V-' 

= :T® 

Note that we have no compensating 5tu-term on the left-hand side of (I30p . To 
circumvent this problem we rewrite the term by using the discrete integration 
by parts formula in time 

_ „k 1 _ uk 1 

'T 'r 


k=l 


k=l 
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Together with the boundedness of c and d (see (A2)) we find: 


= - 


/■*"/■ ^ ^ c(x)Ce(M) - c(x)Ce(u) 

/ / eyn) : -=-dxds 

/o Jo. 


+ f e{v{t)) : c{x{t))Ce{u{t))dx - f e(u°) : c(x°)Ce(u°) dx 
Ju J Q 

= — f f £{v) : -^Ce(u)dxds— f f e(;u) : c(x)Ce(T) dx ds 

Jo Jn Jo Jn ~ 

+ [ e{v{t)) : c{xit))C£(u{t)) dx - [ e(x°) : c(x°)Ce(u°) dx 
Ju J Q 

> -C||c'||i,oo f f |e(;u)||9txl|e(«)|dxds-C' [ [ |e(2;)||e(T)| dxds 

Jo Jn Jo Jn 






C'I|c||l°° [ |e(x(t))||e(u(t))|dx-C'||c||Loo||e(u°)||i2||e(u°)||i4. 
Jn 


_'7^{1)3) 

“•^1 


By using Holder’s and Young’s inequalities, uniform boundedness of \\dtx\\L^{o,T-,L'>‘)i 
Mu)\\l^{ 0 ,t-,l^), l|e(^)llL2(0,T;L2), II^Ct) IIl 2 (o,T;L 2 ) and ||e(T) 11^2(0,T;L2) (see First a pri¬ 
ori estimates), we obtain 


’ < l|5iXllL2(0,tr;L"‘) ll^(T)llL2(0,iT;L"‘) ll^(^)ll-L°°(0,tr;L2) 

< Cs + <^lk(T)|l2,2(oJ^;L4), 


rr' < jlkfe) 


1^2(0,T;L2) + 2 ll^^'^)lli2(0,r;L2) < (J, 


Tr^<csiie(u)iit^(o,T;L-)+mmiih 


The term T| can be estimated as follows: 


rf ^ = 


(d(x)De(u) : e(u) — d(x)De(;u) : £(u)) dxds 


+ / / —d(x)De(u — u) : e(u — u) dxds 

Jo Jn 2r 

f f ^ (d(x)De(^) : £(’«) - d(x)De(i;) : ^(u)) dxds 
Jo Jn 

' -V-' 

+ f j Y “ ^(t)) De(u) : e(i;)dxds 

i/ 0 J o 

'-V-' 

_.7.(2.2) 


+ I / —d(x)De(u — u) : e(u — u) dxds . 
Jo Jn 2r 
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For further estimations we make use of the Ladyzhenskaya’s inequality (see [25]) 

valid for all w € IL^(Q), (32) 

which is a special version of Gagliardo-Nirenberg inequality in 2D (see m)- This 
inequality naturally generalizes to M™'-valued Sobolev functions. 

By using ([32]) . the property d(-) > r/ > 0 and the Lipschitz continuity of d (see (A3)), 
we obtain 


t-rjl 


= E / h (d(x'')De(/) : e{v^) - d(x'=-i)De(/-i) : dx 

= (d(x)De(F(t)) : e{v{t)) - d(x°)De(u°) : e('(;°)) dx 

> : ^iv{t)) - ^d(x°)De(u°) : e(u°)^ dx, 

C [ [ |<9txl|e(T)Pdxds 

Jo Jn 

C [ \\dtx\\LAHv)\\L4^{v)\\L^ds 

Jo 

>-Si f \\dtx\\l4Mv)\\l2 ds - Cs, f \\£{v)\\l4ds 

Jo Jo 

>-Si [ \\dtx\\Hi\\£{v)\\l 2 ds - Cs, [ ||e(T)llL 2 ||e(T)llHi ds 

Jo Jo 

> [ l|5txllii||e(T)|li2ds-C5iC52||e(u)||^2(oy^.i2) -Cl5i(52||e(u)||^2(oy^.Hi)- 


2^(2j2) ^ _ 


> -( 


Note that by choosing <5i = <5 and 62 = and boundedness of l|e(T)lli 2 (oy^;L 2 ) by 

the First a priori estimates, 

rf’2) > ||9ixll^i|k(T)lli2ds-5||e(u)||2,(^-^.^,^ -C,. 

(2 2 ) 

Please notice that in order to treat the term ’ ' in the sequel it will be crucial 
to have established the boundedness of \\dtx\\L^{o,T-,H^) (see First a priori estimates) 
which is due to the higher-order viscosity term —Axt in the damage equation. 

The term r| can be treated by using the Neumann condition (|12cl) and by applying 
the discrete integration by parts formula in time 


j.(3) ^ _ 


f’l 


dtv • 6 dx ds 


/o Jr 

/7 


V ■ dtb dx ds — 


/o Jr 

fir 


J^v{t) ■ b{t) 


dx -b J -b^ dx 


> - 


|ll||L2(r;R")||'9t6||i2(r;R") ds — ||'l^(i)||L2(r.]Rn)||6(t)||2,2(p.Rn) 
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~ ll'f^°llL2(r;R")ll^°llL2(r;R") 

By using the trace theorem L^(r;M”') and the boundedness of 

(see First a priori estimate) as well as of P(^)lli 2 (r;Rr 

||u°||^i and ||^°||^ 2 (r.Kn), we obtain 


t(3) 


- 2 2 C'5||6(t)||^2(r;R^) 

- \h^\\m - ^ll^°lli2(r;R-) 

> -Cs-6\\v{t)\\l,. 


To {T 2 ): With the help of Young’s inequality, we estimate T 2 by 


T 2 > 


^ [ [ d(x)div(De(u)) dxds-C^ f [ c'(x)Vx • Ce(u) 

Jo Jn Jo Jn 


dx ds 




= :T, 


( 1 ) 




-Cs 


10 Jn 


:(x) div(Ce(ri)) dxds+ f f d'(x)Vx • De(u) 

Jo Jn 


dx ds 


= :TP^ 




(33) 


By the following elliptic regularity estimate which follows from [8l Theorem 3.4.1] 
(remember that T is a C^-boundary by (Al)) 

1^11^2 < <^(ll div(De(u;))||2 2 + ||u;||^i + ||De(u;) • z2||^i/2(p.R„)), 

valid for all w G by the Neumann boundary condition (I12cp (remember 

that d(-) > r/ > 0 by (A3)) and by the boundedness of l|T|| 2 , 2 (o_r;/i'i) (see First a priori 
estimate), we obtain 


rpW 


> 


C f f I div(De(x))|^ dxds 
Jo Jn 


^ ^ll^llL2(0,tr;i^2) ^ll^llL2(0,i^;_H'l) C'||De(x) • i^lli2(0,t^;iyl/2(r;Rn)) 


— ^lbllL2(o,t^;J^2) 
= ^\\^\\l^(0 ,U-,H^) 


-c 


^ll^llL2(0,t^;_H-i) 




lL2(0,t.r;_ffl/2(r;R")) 


|C(X) 


d(^)ilL2(0j.;/fl/2(r;R"))W||j^_^ 


+ ||^C£(n) -1211^2(0^ 


tT;iyi/2(r;Rr>)) 


+1 


n(l,l) 


n(l,2) 


The constant C > 0 does depend on rj from (A3). The well-known trace theorem 
yields Ff^/^(r;R”) with a continuous right inverse ^l/2(r.]Rn) 

M”) (see [441 Theorem 8.8]). In the following, we denote the extension of b also 
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by b. We obtain by using the trace theorem, the Gagliardo-Nirenberg type inequality 
the assumptions in (A2) and (A3) and the boundedness of \\^x\\l°°{o,t-,l^) 
\\£iu)\\L°°{o,T-,L^) (see First a priori estimates) 


-‘2 


S + 11^*111^(0,+ lll*ll^^l|ll=(0,!.;L»)) 

S '^(Plll7(o,!,;«.) + r ll*lli» II^IIh- IIVxIIl^IIVxIIh. do) 

^ 0 

“ 2 ~ 2 2 

0 


as well as 




< 


< 


< C(^\\e{u)\\j^2(^Q^^.^ + llk(^)ll^xlll2,2(o^i,^;L2) + ll^(Ce(w))||2,2(o^i,^;L2) 

<^(1 + lk(^)llL4||Vx||^4 ds + \\u\\L^(^o;i.r-,H^)} 

^ lk('“)llL2||e(tt)||_H-il|Vx||2,2(Q)||Vx||j^i ds + ll«lli2(oj,^;_H-2) 

e(^)ll|/i + llVxIll^i) ds + ||«|li2(oj^;H2) 


<C 1 + 


<c 1 + 


/ 

^*"(INI^2 + ||Vx||^i)ds). 


We estimate the remaining terms in (I33p by using again the Gagliardo-Nirenberg 
type inequality (|32]1 . the boundedness of c, d and d (see (A2) and (A3)), and the 
boundedness of ||Vx||j^oo(o,T;L 2 ) and \\£{u)\\l°°{o,T\L'^) (see First a priori estimates) 


Tf <G r(||e(n)||^i + ||Vxf^i)ds, 

Jo 

tP <C r \\u\\]j2 ds, 

Jo 

<c rmwumwhds 

Jo 

f ||Vx||L 2 ||Vx||j^i||e(F)||i 2 ||e(F)||^^i ds 
Jo 
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To (Ts): It can be seen by integration by parts and from the definition of Ip (see 
Definition [231 (hi)-(iv)) that 


^3= r [ ep{dtx)\ydtx\^dxds> 0 . 

Jo Jn 


- To (T 4 ): The term T 4 can be treated by applying the Gagliardo-Nirenberg type in¬ 
equality ([32|) . by using the boundedness of C 4 and of C 2 (see assumption (A 2 )) and 
by using boundedness of ||e(w)||L°°(o,r;L 2 ) (see First a priori estimate). We obtain 


Ti> - 6\\Adtx\\l2(o^t^.L^) - Cs ||e(u)||^4ds 

> -^\\^dtx\\l2^o,t^.L^)-Cs ||e(li)|li 2 ||e(u)||^ids 


- To (Ts): We find by Young’s inequality, by the Lipschitz continuity of f (see (A4)) 
and by the boundedness of ||xIIl 2 (o,T;L 2 ) (see First a priori estimates): 


T 5 > (5||A9tx||^2(oj.^;L2) C'5(||x||^2(o^7’;i2)-I-1) 

— “ '^ll^^iXllL 2 ( 0 j^;L 2 ) — Cs- 


In the following we use the estimates (by the fundamental theorem of calculus and Holder’s 
inequality) 




vP + 

0l|2 


dtu{L)dL 




< lln' 


0||2 


m + ^^ll'^lli2(o,s.r;H2)> 


llx(^)ll|. 2 < 11 x 11 ^ 2 +S.||5*x|li2(0,,,^^,2). 


Now we conclude by taking the above estimates into account; 


(34a) 

(34b) 


Ti > ??||u(t)|||i -T»5-(5(||n||^i(o j^.^ 2 ) + ||u(t)|||i) - 5 ||9tx(s + t)||^i ||u(s)||^i ds, 

T 2 > ^lkllHi{ 0 ,i.r;i^ 2 ) -Ds-Cs f (^11x11^2 + 11^11^1(0, s^;/f 2 ) +G||T ||^1 1194x11^2(0,ds 

^ 0 

P^T _ 

— '^II'“II_H- 1 ( 0 , 4 ^;_H' 2 ) ~ ^ ll^llHl/ 2 (r;R") IIxII_H -2 

73 > 0 , 

^4 > ~ ~ ll'“'llHi(0,Sr;H2) ds - (5||xll^l(0,t^;H2), 

> - Cs - <5||xIIj^1(o,4.^;H2)’ 

where the constant > 0 continuously depends on (besides 6) 

C5 = C4||u°|| HX lb°llL2, IIx^IIj/X ll^llL2(o,T;I,2(r;R")), PIIl2(o,T;L2)) 
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and the constant Ds > 0 continuously depends on (besides 6 ) 


Ds = Ds{\\u'^\\ ||u°||jyl, ||x°IIh 2) ll^llL2(-o^T;Hl/2(r;R"))ni^l(0.'r;i^(r;R"))’ PII-L2(0,T;L2))- 

By adding the identities (IHU]) and (f 3 T]l . using the estimates for Ti,..., T5 developed above, 
using the i 7 ^-regularity estimate (see [H Theorem 3 . 4 . 1 ]) 

11^11^2 < C'(ll^'»^lli2 + \\w\\hi) valid for all w € 

applied to x dtx (note the boundary conditions in ([HD), we obtain 

\\v{i)\\Hi + \\W)\\h2 + + llxll|-i(0,tr;H2) 

+ II div(c(x)Ce(u) + d(x)De(T)) |li2(oj^.i2) 

< Ds + Cs (||xIIh2 + \W\\m{o,sr-,H^) + 

+ '^(11^11^1(0, + ll^(^)llii + II div(c(x)Ce(u) + d(x)De(u))) + ||xll^i(o,i,.jy 2 ) 

+ + PllHi/2(r;R")llxlli2ds. ( 35 ) 


By choosing <5 > 0 small, the first (5-term on the right-hand side of (I35D can be absorbed 
by the left-hand side. Furthermore, for later estimates, 5 should also satisfy 


(5 < 


8(||<9tX||^2(o,T;Hl) 


1 


(36) 


Indeed, the denominator of the right-hand side is bounded from above by the First a priori 
estimates, hence the right-hand side is bounded from below and (5 > 0 can be chosen such 
that ([36D holds. 

We infer from the estimates ([SSP and (I5SD 


k 

ot^ < Ds -\- ^ 
i=i 


(37) 


with 


:= 


7^: 


— Cs -\- 


II (X" 


8(ll^iXll/^2(o,r;Hl) ll^llL2(o,r;Hl/2(r;R")) 


In the following, we will choose a time to > 0 such that for all small r > 0 and all 
k = l,..., (to)^/r: 

0<r7"<^. (38) 
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Indeed, we know by the First a priori estimate that 

M 

< C uniformly in r, 

k=l 

where C > 0 denotes the constant C in (l23]i . Thus 

r||f^||^i < C uniformly in r and in k. (39) 


By choosing 


to : — 


1 

4W’ 


(40) 


we get for all A; = 1,, {to)^/T: 


+ tC5 X ^7^ X llf^ll^i 
<{to).r^*0+'r 

< tCs + o + TCs{to + 'r)||n^||^i 

O 

< tCs g 4 ’^^6^ 

< tCs{1 g + 4’ 

Consequently, for small r > 0, estimate (l38|) is fulfilled. 
Finally, by ensuring (I38p . (|37l) rewrites in the desired form 


and, therefore, 


< 


Ds 


k-l 


^ ^ \ — 7"X' 


'^J 


-a-' 


a <-h > T— w. 

“2^2 

i=i 


We are now in a position to apply the discrete version of Gronwall’s lemma in the sum 
form (see, e.g., m page 26]) and obtain 


a 


< 


Dx 

_ = l ‘ 2 


We obtain boundedness of uniformly in r and k = 1,..., (to)r/E. Therefore, (ii) 
is shown except the boundedness for \\ii3{dtXT,i3)\\L'^(o,(to)T\L'^)' latter follows by a 

comparison argument in (I12bh . 
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To (ii) - global-in-time estimate: 


The main observation to obtain global-in-time estimates is that the local estimates above 
can not only be performed on the time interval [0, (tQ)r] but also, with minor modifications, 
to each interval ^ [0,T] such that {It — < to, where to > 0 from (^01) depends 

on quantities which can be bounded globally in time by the First a priori estimates. Thus 
we find a to > 0 such that the Second a priori estimates can be performed on each interval 
interval [s^,tr] C [0,T] with |t,- — < to- 

To conclude the proof, let 


:= ~ I i € N such that k-^ > rj| 

ir ■■= = min |jr | j € N such that to < Tj|. 


We define the time intervals 


:= [t^,t^ + h]n[o,r] 


for all A: = 0,..., with N := |'T/(to/2)] — 1 where f-] denotes the ceiling function. 

We apply the local-in-time estimates above to each interval and obtain constants 
Co,..., Cat > 0 which continuously depend on 


Ck — Ck{\\Ur,j3{ir)\\H2, ||Xr„3(tr) , 

ll^llL2(o,T;//l/2(r;]R'i))nHl(0,T;L2(r;R"))) II^IIl2(0,T;L2)) , 


k = 


such that for all r, /3 > 0 and all A: = 1,..., 


\'^T,l3\\L°°(Ik-H'2) < Ck, 


\\XT,l3\\H^{Ik-H^) < Ck, 

< Ck, 

\\Xt,I3\\l°°{I^-,H^) < Ck, 
Ul3{dtXT,/3)\\L2{lk-L^) < Ck- 


To obtain a global bound, we can argue by induction. We sketch the argument: 

Suppose we have given the a priori bound Ck-i for the time interval . By definition, 
we hud Thus 

ll'*^T,/3(tr) — IIW,/3||£,oo('jfe-l.JA2) — Ck — l, 

||'yT,/3(tr) Ilii'l — II^^t,/3|Iaoo(7*-1;JA1) — ^k-l, 

II AT,/j(^r) IIh^ — II At,/3 II— Ck—l- 

Consequently, we find an a priori bound Ck > Ck for the solutions on the interval by 

Ck ■= 2/C) l|ft|lL2(o,T;fl'l/2(r;R"))n//l{0,'r;i^(r;K"))> 11^11^2(0,r;L2)). 

1^1 n?/n l-^l —1 
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Note that Ck does only depend on 


ll^llL2(0,T;//l/2(r;K"))n//l(0,T;L2(r;K"))) II^IIl2(0,T;I,2)) . 


□ 

We perform the limit passage r 0 and /? i 0 separately in order to show existence of strong 
solutions for both cases: namely for /3 > 0 and /3 = 0 in Definition 12.51 (i) and (i). The a priori 
estimates give rise to the following convergence properties along a suitably chosen subsequence. 


Lemma 2.10 (Convergence properties) 

There exist limit functions for every (3 > 0 (we will also write u := uq, X ■= Xo) 

U0 G H\0, T; M^)) n T; H\n; M”)) n H^{0, T; L2(D; M^)), 

with 


Uj 3 { 0 ) = a.e. in D, a.e. in D, X/3(0) = in D, 

Uj 3 = b a.e. on S 

such that 


(i) for fixed /3 > 0 and r J, 0 (along a subsequence): 



weakly in F^(0, T;M”)) 
weakly-star in VT^’°°(0,T; 

(41a) 

(41b) 


weakly-star in L°°(0,T;M"')), 

(41c) 


strongly in H^{0,T] 

(41d) 

'^13 

strongly in L°°(0, T; M"')), 

(41e) 

^r,/3; ^r,/3 ^ ^(3 

a.e. in kl X (0, T), 

(41f) 

Vt,p dtup 

weakly in H\0,T; L‘^{n-,R^)), 

(41g) 

XT,f} Xh 

weakly in H^{0,T;H^{XI)), 

(41h) 

^T./3Wr,/3 

weakly-star in L°°{0,T; 

(41i) 


strongly in L^{0,T-, for all 

(41j) 


uniformly on XI x [0,T], 

(41k) 

f,l3idtXT,i3) ipidiXp) 

weakly in L^(0,T; L^(D)), 

(411) 

4- 0 (along a subsequence): 

Up ^ u weakly 

in H^{0,T-,L‘^{Q-,R^))nH^{0,T;H^{Q-,R'^)) 

(42a) 

weakly-star in T; M”)), 

(42b) 
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Up ^ u 

strongly in Lr^(0, T; M”)), 

(42c) 

up ^ u 

a.e. in Q X (0,T), 

(42d) 

Xp^X 

weakly in T ; 

(42e) 

Xp^X 

strongly in L^{0,T; for all /i > 1, 

(42f) 

Xp^X 

uniformly on Q x [0,T], 

(42g) 

CpidtXp) C 

weakly in forage L‘^{0,T, L‘^{n)) 

with ^ G 9/(_oo,o] (^tx) in Q, X (0,T). 

(42h) 


Proof. 


To (i): Properties (|41ap - (l41cp and (l41fl) -( [lTj] ) can be obtained by standard compact embed¬ 
dings, whereas (I41dp . ()41ep and ()41kh can be obtained by the Aubin-Lions type compact¬ 
ness result in [12] (please note that = dtUr^j^)- It remains to show (j41ip . 

By Lemma [2^ we find a cluster point rjp € T^(0, T; L^(f2)) such that along a subsequence 
r 4 , 0 


^p{dtXT,p) Vp weakly in 


(43) 


We have to show = ^pidtxp) to finish the proof. This will be achieved by exploiting 
maximal monotonicity of the graph (also known as Minty’s trick). At first we ob¬ 
serve that due to the monotonicity of viewed as an operator from L^(0, T; L^(Q)) to 
L2(0,r;L2(O))' one has 

{(p{dtXT,p) - ip{v), dtXT,p - > 0 

for all functions v G L^(0, T; L^(n)). In order to pass to the limit r 0 it remains to show 

{^f}{dtXT,p),dtXT,p)< {npidtXp)L2(^L2y (44) 


Testing equation (I12bl) with dtXr^p and integrating over fl x (0, T) in space and time 
and passing r 0 by using weak lower-semicontinuity properties for the ff \dtXT,p\‘^- and 
ff I V5tXT-^y3p-term, we obtain 


lim inf / / 

’■fo Jo Jn 

f- 

> lim inf / 
"^io Jo 


-^p{dtXT,p)dtXT,p dx ds 


/ |<9tXr,/3pdxds-hliminf / / \X/dtXT,p\‘^ dx ds 

Jn ■rfo Jo Jn 


+ lim / VXr,B ■ '^dtXT,pdxds 
Jo Jn 


+ 1™ ^ X Ce(Mr,/3) : 4ur,g) + f'iXr,p)^ dtXT,p dx ds 

> [ [ (\dtXp\‘^+ Ndtxp\‘^) dxds+ [ [ Vxp -"^dtXpdxds 
Jo Jn ^ ' Jo Jn 

^ I In ■ ^iup) + f{xp)')dtxpdxds 
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^ L In ~ ~ ^c'(x/ 3 )Ce(n/ 3 ) : e{uf}) + /'(x/?)^ dtXii dx ds. (45) 

Note that we also get 

-Vp = dtXp - ^Xi3 - A5tX/3 + ^c'(x/ 3 )C£(u^) : e(u^) + f'{xp) a.e. in x (0, T) 

by performing a limit passage r 4- 0 in (|12bp after testing with a function, integrating 
and using the already known convergence properties (j41al) - (l41kl) and (|43D . In combination 
with (145 p and multiplying by —1, we find the desired limsup-estimate (14411 . Thus 

{vp - ^isiv),dtxi} - v)> 0 (46) 

for all functions v G L^(0, T; L^(n)). By maximal monotonicity of = dip viewed as an 
operator L^(0, T; L^(n)) —)> L^(0, T; L^(fl))', we find r/p = Cpidtxp)- 

To (ii): Since the a priori estimates in Lemma 12.91 are also independent of /3, we obtain an 
analogous result for the limit case /3 0. 

It remains to establish (I42hl) . Convexity of Ip viewed as a functional T^(0, T; L^(II)) ^ M 
implies 


IpidtXfs) + {Ci3idtxp),v - dtxn) l2(l2) ^ ^/ 9 (^) (47) 

for all functions v € L^(0,T;L^(fl)). With the same arguments as in (i) we obtain the 
limsup-estimate 


limSUp(^;3(5tX/3),5tX/3)£,2(i2) < {v,dtx)L2(^L2y 

where r/ is a cluster point of {5tX/3}/3G(o,i) ™ weak topology of L^(0, T; L^(n)). 

Now let V € L^(0, T; L^(n)) with u < 0 a.e. in II x (0, T) be arbitrary. By using Ip{v) = 0 
and Ip{dtXi3) > 0, we obtain from (I47p after passing to /3 4- 0 for a subsequence 

(r/,u - 9tx)^2(i2) < 0. (48) 

We also infer from (I47p the following boundedness with respect to /3: 

If^idtXp) < C. 


For the limit /3 4- 0 we thus obtain 

dtX < 0 a.e. in 14 x (0,T). (49) 

Estimate (flSP together with ()4UP yields r] € dI^ 2 (^Q rp.i 2 (Q))(9tx) with the indicator function 
/^ 2 (q 2 ’;l 2 : L^(0,T; L‘^{Q)) —>■ M U {+oo} and := {v G | u < 0 a.e.}. 

We end up with rj G dI(^_aofi]{dtx) a.e. in 14 x (0,T). □ 
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Theorem 2.11 Let the Assumptions (Al)-(A4) be satisfied and the data ,b,€) from 

m and (USD be given. The following statements are true: 

(i) Regularized case (fi > Oj; There exists a strong global-in-time solution {up,xp) lb,e 
sense of Definition \2. 5\ (ii) which satisfies (fTH|l . 

(a) Limit case (fi = Q): There exists a strong global-in-time solution {u,x) in the sense of 
Definition 1 2. 5\ fi). which satisfies (fTHl) . 

Proof. 

To (i): By multiplying the systems ()12ap . ()12bp . ()12cll and (I12dp with test-functions, inte¬ 
grating over space and time, we may pass to the limit r 0 for fixed /? > 0 by utilizing 
Lemma 12.81 and Lemma 12.101 (i) and standard convergence arguments. Then, switching 
back to an {x,t)-a.e. formulation, we obtain a strong solution of system (fT P -(|5 p . 

To (ii): The transition /3 4- 0 can be conducted as in (i) by utilizing Lemma 12.101 (ii). □ 


Remarks to the proof of Theorem 12.111 

(i) The regularizing term —Axt in (I12bl) is needed in order to obtain an {H^)-honnd for 
X in the first estimate. This, in turn, was particularly necessary to estimate 6 in (I.SOp and 
to estimate the term T/^ in the second estimate. 

(ii) In the mathematical literature the elasticity equations (ll2aD is sometimes tested with the 
function — div(De(tti)) to gain higher-order estimates for u (see [5U1133 SO] )• However, 
due to the nonhomogeneous Neumann boundary condition (l8cP in our case, it is more 
convenient to test with — div(c(x)Ce {u) + d(x)De(nt)) since, otherwise, integration by 
parts in space of the term — ff uu ■ div(De(Mt)) yields unpleasant terms even after using 
the boundary condition (l3^ (cf. estimates for Ti in the second estimate). 

2.4 Continuous dependence on the data 

We are going to show continuous dependence on the data of the strong solutions of the PDE 

system given in Definition 12.51 fil and (ii). 

Theorem 2.12 (Continuous dependence) Let the Assumptions (Al)-(A 4 ) be satisfied. More¬ 
over, assume that d = 1 in (A 3 ) and one of the following condition: 

• Let (ni, xi) and {u2, X2) be both strong solutions according to Definition \ 2 . 5 \ (ii) (for fi > Oj 
with data {ui,Vi,Xi,bi,ii) and (n®, n®, X 2 ) ^ 2 , ^’ 2 )- 

• Let (ni,xi) o,nd (02^X2) be both .strong solutions according to Definition \ 2 . 5 \ (i) (for fi = d) 
with data (n)*, n)*, x?, &i,^1) and (^2, ^2, X2) ^2,^2)- 

Then, 

ll^^l — '*^2||wC°°(0,T;L2)n//l(0,T;//l) + IIXl “ X2H^l(0,T;Hl) 

< C^lln? — + 11^1 — + llxi — + IKl “ ■^2 ||l2(o,T;L 2) -I- ll&l — f>2||L2(0,T;L2(r;K"))^ , 

(50) 
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where the constant C > 0 continuously depends on 

C = c[\\ui\\u, \\u 2 \\u, IIXiIIa’, ||X 2 ||A’y 

Proof. For notational convenience, define 


U\=Ul- U2, 

£ := £i — £ 2 - 


X-=Xi- X2, 
:= X? - X2> 


:= Ui — U 2 , 

b:=bi- 62 , 


Let t € [0, T] be arbitrary. Firstly, testing the damage equation (jlbli for each solution with Xti 
subtracting the resulting equations and integrating over XI x (0, t), we obtain 

[ j (lx/ + Vx • Vxt + |Vx/ + ^c'(xi)(Ce(ui) : e(ui) - Ce(u2) : e{u2))xt) dxds 

+ / / (^(c'(xi)-c'(x2))Ce(u2) :e(n2)Xi + (/'(xi)-/'(X2))Xt + (6-6)Xt)d2:ds 

o 0 t/ ri 


= 0 . 

By assumption, we know 

€ 9/(_oo,o](9tXi)> * = 1,2 
Ci = I'pidtXi), i = l,2 

It follows from the monotonicity of 9/(_oo,o] ^-^^d (see Dehnition l2.3jl . respectively, that 

(Cl - 6)%* = (Cl - C2)(xi - X2)i > 0. 

Therefore, by 


(51) 


if /3 = 0, 
if /3 > 0. 


\\Xt\\l^o,t;m) + ^l|Vx(t)|li 2 - ^||Vx 1 li 2 


0||2 


< - 


1 / 


c (xi)Ce(ri) : e{ui + U2)Xt d^c ds 


0 Jn 


= :Ti 





o(c'(xi) - c'(x2))Ce(u2) : e(u2)Xidxds 
lo Jn ^ 

V 

= .T2 


U'ixi) - /'(X2))Xtdxds. 


(52) 


0 Jn 


=.Tz 


By using Holder’s and Young’s inequalities as well as standard Sobolev embeddings, the Lipschitz 
continuity of d (see (A2)) and the Lipschitz continuity of /' (see (A3), we hnd 


Ti < 


C'||e(ui+U2)||Loc(o,i;z,4) [ ||e(w)||L2|IXtllL4ds 

Jo 
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T 2 < C||c'||Lip||e(7X2)|lioc(o,t;L4)^ \\x\\L4Xt\\L'tds 

< 4Xt\\L2{0,t-,m) + ^511x11^2(0, 

73 < ||/'||l,ip / ||xllL2||XtllL2ds 

Jo 

< '^llXt|li 2 ( 0 ,t;L 2 ) + C'<5||x||i2(0,t;L2). 

Applying the estimates for Ti, T2 and T3 to ([ 5 ^ . we obtain 

\\Xt\\L^{o,f,m) + I|Vx(^)IIl 2 < C' 5 (^||x°||^i + 11x11^2(0,t;Hi) + \\4\h{o,t-,m)) + 4Xt\\h{o,t-,m)- 

(53) 


Secondly, we test each of the corresponding elasticity equations (flail for ui and U2 with Ut 
and obtain by subtraction and integration over x (0,t): 

J j [uttut + c(xi)Ce(u) : £{ut) + (c(xi) - c(x2))Ce(n2) : £{ut) + Defri*) : £{ut)^ dxds 



0 JV 


&'Utdxds+ / / i-Utdxds. 

Jo Jn 


This implies 

\\\ut{t)\\l2 - ^Il^°lli 2 +r/l|e(^^t)lli 2 ( 0 ,t;L 2 ) 


< - 


n c(xi)Ce(ri) : e(iit)dxds- f f {c{xi) - c{x2))C£{u2) : £{ut) dx ds 
! Jo Jn 


= :T4 


= -T5 


+ f j h-utdxds+ r [ 
Jo Jr Jo Jn 

= :n 


£ ■ Ufdx ds 


(54) 


0 Jn 


= :T7 


Standard estimates yield 

T4 < 6\\e{ut)\\‘j^2(^Q + Cs\\£{u)\\‘j^2(^Q 

^5 < C'llxllL2(oy;L3) ||e(tt2)||L°°(L6) ||e(rit) 112,2(0,f;L2) 

< + ^5||xlli2(o,i;22i), 

Tg < <^||'lit|li2(0,t;L2(p.Rn)) + C'(5||^>|li2(0,i;2,2(p.]Rnp 

< 4'^t\\\,^{0,t\H^) + ^ 5 ||h||p 2 ( 0 ,i; 2 , 2 (p.]Rn)), 

T7 < <^||'W't|li2(0,t;L2) + C'<5||^llL2(0,t;L2). 

Applying the estimates T4, T5, Tg and T7 to (IMD shows 

ll'“•^(^)llL 2 + lk(^i)llL 2 ( 0 ,t;L 2 ) 
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( 55 ) 


< C5{\\v^\\\2 + ||w|li2(o,j;Hi) + ll^llL2(0,t;Hi) + ll^llL2(o,t;I,2)) 

Adding (fM|l and ([55]), we obtain 

lk('“t)llL2(o,4;L2) + ||Atlli2(o,t ;Hl) + \\Ut{t)\\l2 + \\Vxml2 

< Cs(\\v°\\‘i2 + ||A°II^1 + IIaIIl 2(0,J;H1) + Il'“lli2(0,t;//1) + ll^llL2(0,t;L2) + 11 ^ 11 ^2 (Q^t ;L2 (p-Rn)) ^ 
++ ll^tllL 2 ( 0 ,t;Hl)) • (^ 6 ) 

Now, adding ||w||^ 2 (o,t;_H'i) + \\'^t\\L^(o,t-,L^) using 

llAWIli2 = IIa°+ f Xt(s)ds||i2 < C(||x°|||2 + ||Xt||i2(0,t;L2)). 

^ 0 

and Korn’s inequality 

Iklki < C(||u ;||^2 + ||e(?u)||i 2 ) 

holding for all w G H^{Q,]W^), the estimate (|56l) becomes 

+ IIXtllL2(0,t ■,m) + Ik^wiii2 + iix(t)ii^i 
< C'<5(ll'*^°lli2 + IIx°IIhi + ll^lli2{o,4;Hi) + ll'“llHi(o,t;L2) + \\'^\\h{o,t-,m)) 

+ C'5(^||-^||i2(o,t;L2) + ll^lli 2 (o^t;i 2 (r;R»»))) + + IIXtllL 2 ( 0 ,t;Hi)) ' 

By choosing <5 > 0 small and noticing 

Il'“lli 2 ( 0 ,t;//l) <c(^\\u°\\jji + Il'“t|li2(0,s;_f/l) 

we get 


\\'^\\‘H^iO,t;H^) + ll^tllL 2 ( 0 ,t;Hl) + ll'“t(^)llL 2 + IIx(^)IIh 1 

< c(\\u^\\j{l + \\v°\\l2 + ||A°||^i + ll^llL2(0,t;L2) + ll^llL2(0,t;L2(r;Rn))^ 

lo Il'“illi 2 ( 0 ,s;// 1 ) + llxll^l) ds. 

The claim follows by Gronwall’s lemma. □ 

The continuous dependence result in Theorem 12.121 as well as the a priori estimates in Lemma 
12.91 yield the following corollaries. For notational convenience we will make use of the spaces U 
and A defined in Subsection O (see m)- 

Corollary 2.13 (Uniqueness) Strong solutions in the sense of Defi,nition \2.,’^ (i) or (ii) with 
constant viscosity D are unique to given initial-boundary data {u^,v^,x^,b,i). 
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Corollary 2.14 (A priori estimates) A strong solution {u,x) € W x A for the system in 
Definition \2. 5\ (i) or (ii) with constant viscosity B and given data satisfies the a 

priori estimates 

\W\\u < C, IIxIIa- < C, II?IIl2(Ox(o,t)) < c, 

where the constant C > 0 continuously depends on 

C = C[\\u^\\h2, ||v°||_f/i ) ll^llL2(o,T;_H'l/2(r;R’^))nHl(0,T;L2(r;R"))) II^IIl2(o,T;L2)) • 

3 Optimal control problem 

In this section we establish the announced optimal control problem for the damage-elasticity 
system m-m- From now on, we assume for the viscosity tensor B = D, i.e. d = 1, in order to 
apply the well-posedness result from the last section. 

Let U, A and B be given as in Section 12.11 Our aim is to approximate with y prescribed 
damage profiles by controlling the Neumann boundary data b £ B for the stress tensor a in 
(j3bl) . The cost functionals measures the deviation from the prescribed profiles at the final time 
or/and at all times during the evolution in an L°°-norm. We make the following assumptions: 

(Ol) We assume that Xq, Aq and As are given non-negative constants which do not all vanish. 

(02) The target damage profiles are given by 

XQ E L~(Q), XT G 

(03) The admissible set of controls Badm ^ B is assumed to be non-empty, closed and bounded. 
Remember that B := L^(0, T; M”)) n 1/^(0, T; L^(r; M”)) (see Subsection 12. ip . 

Remark 3.1 A typical choice for Badm would be 

Badm — G B I b^in b ^ b^aiax a.c. in B and ||f ^||/5 ^ ALj-, 

where M E (0,oo) denotes the maximal B-cost and bm,in,bmax G B the minimal and maximal 
cost functions satisfying bmin < bmax o-.e. in S. 

We define the following tracking type objective functional 

J{Xib) := -^\\x — Xq\\l^{q) + -i^WxiT) — Xt\\l°°{q) + ^II^IIl2(S;R' 
where our overall optimization problem reads as 
minimize J (x, b) over A x Badm 

s.t. the PDF system in Definition 12.51 (i) is satisfied for an u £U. 


(57) 

(Po) 


Remark 3.2 Let us emphasize that we may also choose || • ||| 2 -^erms instead of the || • \\l°°- terms 
in the cost functional ill). The existence results presented in this section work for both cases. 
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We recall that the system m-m is an initial-boundary value problem, which admits by Theorem 
12.111 and Corollary 12.131 for every satisfying (fT7D and (fTS|l a unique solution 

(n, x) G Z// X T’ in the sense of Definition 12.51 (i). Hence, the solution operator 

^ {u,x) 


with 


1 := b) I satisfying (fT7|) and ifTH]) }. 

0:=UxX 


is well-defined. Moreover, for fixed data {u^,v^,x^,^) the control-to-state operator 

'S'o : Sadm O, 6 I—> (u, x) 

is also well-defined, and the optimal control problem ( |PoD is equivalent to minimizing the reduced 
cost functional 

m ■.= J{So\2{h),h) 

over Badm, where S'o |2 denotes the second component of Sq, i.e. Sq = (Poll)'S'o| 2 )- 

For /3 G (0,1), let us denote by Sjs the operator mapping the control b G Badm into the 
unique solution (u/ 3 ,X/ 3 ) G O to the /3-regularized problem in Definition 12.51 (ii). 

Remark 3.3 In view of the continuous dependence result in Theorem \2.1^ the operators Tq 
and So are well-posed in a larger target space (VF^’°°(0,T; L^(D;M"')) n P^(0,T; P^(D;M"'))) x 
whieh eontains the spaee O. This fact is important for the sensitivity analysis 
of these operators. But in this section, we are interested only in existence of optimal controls, so 
this result is not needed. The sensitivity analysis which also establishes the optimality conditions 
of first-order will be treated in a forthcoming paper. 


3.1 Existence of optimal controls to ([^| ) via /?-regnlarization 


The following lemma is the basis for the main result in this section. 


Lemma 3.4 We have the following continuity properties: 
For a given sequenee {fe/ 3 }y 3 g(o,i) C B and b G B with 


bp ^ b weakly in B as /3 4- 0, 


it holds 


Spibp) ^ So(b) 

Sr,ibp) -G Sr,(b) 

Soibg) ^ Soib) 


weakly-star in O as /3 4- 0, 

weakly-star in O as /3 fO for every p G (0,1), 

weakly-star in O as /3 4- 0. 


(58) 


(59a) 

(59b) 

(59c) 


Proof. Let {up,xii) = Sp{bp). Then, {up,xii) is a solution to the /3-regularized system in the 
sense of Definition 12.51 fiil with Neumann data bp. Since {bp} C H is bounded by ([581) . we obtain 
the a priori estimates from Corollary 12.141 In particular. 


Sp{bp) —>■ (tt,x) weakly-star in O as /3 f 0 (for a subsequence). 
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for some {u, x) ^ 

We see from the proof of Lemma I2.1UI that the convergence properties in Lemma 12.101 (ii) 
hold for a subsequence /3 | 0. By using these convergence properties as well as we can 
pass to the limit for a subsequence in the /3-regularized PDE system (lllal) - ()lldp (cf. the proof 
of Theorem 12.lip . We obtain that {u,x) satisfies the limit system in Definition 12.51 (i) to the 
Neumann data b. In other words, So{b) = {u,x)- By the uniqueness of solutions shown in 
Corollary 12.131 we see that Sp{bi 3 ) convergences weakly-star to (u, y) for the whole sequence 
/? I 0. Hence, (|59al] is shown and (I59bp and ()59cl) follow with the same reasoning. □ 


Corollary 3.5 We also have the property 

limj(5^|2(6),6) = J(5 o|2(6),6) 

P-l-U 

for allbG Badm- 

Proof. Set bp := b and apply Lemma IHITI In particular, (I59al) implies 

5'/3|2(^) ‘S'o| 2 (&) weakly-star in T as /3 4- 0- 

Define Xh ■= Sp\ 2 {b) and y := 5'o|2(f))- A standard compactness result shows (see [32]) 

Xfi^X strongly in C°{Q) as /3 3 0 

and the claim follows. □ 


Theorem 3.6 Suppose that Assumptions (Al)-(A4) as well as (01)-(03) are satisfied. Then 
the optimal control problem dPoP admits a solution. 


Before proving Theorem 13.61 we introduce a family of auxilliary optimal control problems (P/?), 
which are parametrized by /3 G (0,1). We define 


minimize JT"(x, b) over X x Badm 

s.t. the /3-regularized PDE system in Definition 12.51 (iil is satisfied. 


(Pp) 


The following result guarantees the existence of an optimal control to {Pp). 

Lemma 3.7 Suppose that the Assumptions (Al)-(A4) as well as (01)-(03) are fulfilled. Let 


j3 > be given. Then the optimal control problem (Pp) admits a solution. 


Proof. Let {6"'}nGN C Badm be a minimizing sequence for (Pp), and let {u'p,x'p) = SpifP), 
n € N. By the boundedness and closedness of Badm (see (03)), we find a function b G Badm and 
a subsequence of {6”'} (we omit the subscript) such that 


b^ ^ b weakly in H as n 3 oo. 


Lemma 13.41 yields 


Sp{b^) Sp{b) weakly-star in O as n t oo. 


35 
























Let {ui3,X/3) ■— S/sib)- We particularly find 

= Sf^^ 2 {b^) 'S'/ 3 | 2 (^) = X /3 weakly-star in as n t oo. 

A standard compact result reveals 

strongly in as n t oo. 

It follows from the sequentially weak lower semicontinuity of the cost functional J that h is an 
optimal control for (Pp), i.e. 

J{X0,bp) 


□ 


Proof of Theorem 13.61 By virtue of Lemma 13.71 for any /3 € (0,1), we may pick an op¬ 
timality pair 

(X/3) bp') € A X 13adm 


for the optimal control problem \Pph- Obviously, we have {up,xp) = Sp{bp), /? E (0,1). 
assumption (03) and Lemma [3~^ we find functions (u, x) € O and b E Badm with So(b) 
such that 


By the 

= {u,x) 


iup,Xp) (u,x) 
bp ^b 


weakly-star in O, 
weakly in B 


(60a) 

(60b) 


as /3 0 (for a subsequence). 

It remains to show that (Xj &) is in fact an optimality pair of ( |PoD . To this end, let b E 
Padm be arbitrary. In view of the convergence properties (|6n|l and the sequentially weak lower 
semicontinuity of the cost functional, we have 


J{x,b) < limmij{xp,bp) = limini J{Spi 2 {bp), bp). 


By using the optimality property of (Pp), we obtain 

limi)nf.7(5’^|2(&/3),&/3) < liminf J"(5^|2(6), 6). 
Finally, the convergence property in Corollary 13.51 shows 


hininf J(5^|2(6),6) = J(5 o|2(6),6). 

In conclusion, we have proven i7(5'o|2(^),^) < {So\ 2 {b), b)■ □ 

Remark 3.8 Theorem I,?, hi can also be shown in the spirit of Lemma \3.1\ However, the proof 
presented via convergence of 13-approximations might be of interest in view of the implementation 
of optimality systems. 
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3.2 An adapted optimal control problem to (j^]) 


Theorem 13.61 does not yield any information on whether every solution to the optimal control 
problem ( |Po[ ) can be approximated by a sequence of solutions to the problem {P/s). As already 


announced in the introduction, we are not able to prove such a general “global” result because 
of the lack of uniqueness of optimizers. Hence we cannot guarantee that the weak limit of the 
sequence of optimizers to problem (P/s ) converges always to the same optimizer of ( IPqI ). Instead, 


we can only give an answer for every individual optimizer of ( |PoD , which is the reason why this 
solution is called a “local” result (see also Remark 13.111) . For this purpose, we employ a trick 
due to [21 Section 5 - Proof of Theorem 1]. 

To this end, let ((u, x),6) € Ox Badm-, where (u, x) = So(b), be an arbitrary but fixed 


solution to ( Pq ) • We associate with this solution the adapted cost functional 


1 , 


J{Xi b) := J{Xi b) + -\\b — b\\L2/Y:;Rn.) 
and the corresponding adapted optimal control problem 
minimize J (x, b) over X x Badm 

s.t. the /^-regularized PDF system in Definition 12.51 (hi is satisfied. 


{Pfi) 


With a proof that resembles that of Lemma 13.71 and needs no repetition here, we can show the 
following result: 

Lemma 3.9 Suppose that the Assumptions as well as (01)-(03) are fulfilled. Let 

/3 € (0,1) he given. Then, the optimal control problem (P/s) admits a solution. 


We are now in the position to give a partial answer to the question raised above. More precisely, 
we show the following theorem: 

Theorem 3.10 Let the Assumptions (Al)-(A4) and (01)-(03) he satisfied. Suppose that {x,b) G 
31 X Badm is any fixed solution to the optimal co ntro l problem dPoD . Then, there exists a pair 
iX/ 3 ^b/s) £ X X Badm Solving the adapted problem (P/s) such that JifX/sJ^h) J^ix^b) as j3 fQ. 


Proof. For every (3 G (0,1) we pick an optimal pair (x/s^b/s) G X x Badm for the adapted 
problem (P/s). By the boundedness and closedness of Badm (see (03)), there exists a b £ Badm 
satisfying 


b/s ^ b weakly in P as /3 4- 0. 

Owing to Lemma 13.41 we find 

{ujS.Xp) = S/s{b/s) So{b) =: {u, x) weakly-star in O as /3 I 0. 
and, particularly, 


(61) 


X/s^X strongly in C'^{Q) as /3 | 0. (62) 

We now aim to prove that b = b. Once this is shown, we can infer from the unique solvability 
of the state system (see Theorem 12.12p that also (tt,x) = iu,x)- 
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Indeed, we have, owing to (l6T|) . ([62]), the sequentially weak lower semicontinuity of J, and 
the optimality property of (x, &) for problem (|PoD, 


1 , 




/34.0 


>t7’(X)&)+ „\\b ^llL2(5].Rn) 


(63) 


On the other hand, the optimality property of (x/3, bjs) for problem (P^) yields that 

= JiSfsi2ih)^h) < JiSf,^2{b),b)- 

Whence, taking the limes superior as /3 4- 0 on both sides and invoking Corollary 13.51 we find 

limsup J(x/3,&/3) < J{So\ 2 (b),b) = Jix,^ = J{x,b). (64) 

/ 3 i 0 


We obtain by combining (l63l) and ([Ml) 

2 11^ “ ^IIl2(S;R'i) = 0- 

Thus b = b and, consequently, (n, x) = (u, x) by Theorem 12.121 

Finally, by using 6 = 6 in (l63l) and ([6^ . we end up with lim^4,o <JiXp:bp) = J{x^b). □ 


Remark 3.11 As we have seen in Theorem. \S.fk every weakly convergent subsequence of mini- 
mizers of {Pp) converges to a minimizer of ( |PoD . However, since the problem dPpP might not 

might converge to different 


be uniquely solvable different subsequences of minimizers of 
minimizers of (|PoD. By considering the adapted control problem (P/ 3 ), we force the minimizers 


of (P/3) to converge to the desired minimizer of ( |PoD as shown in Theorem \d.l(A For further 
details to non-convex optimal control problems we refer to 
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